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Lattice QCD determinations appropriate to hadron spectroscopy are reviewed with emphasis on the glueball 
and hybrid meson states in the quenched approximation. Hybrids are discussed for heavy and for light quarks. 
The effects of sea quarks (unquenching) are explored. 



1. INTRODUCTION 

Quantum Chromodynamics is generally ac- 
knowledged to be the theory of hadronic inter- 
actions. Its perturbative features are well un- 
derstood and provided the main motivation for 
its adoption. To give a complete description of 
hadronic physics, it is essential to develop non- 
perturbative methods to handle QCD. QCD is 
a quantum field theory which needs to be reg- 
ulated in order to have a well defined mathe- 
matical approach. Inevitably such a regularisa- 
tion will destroy some of the symmetries of QCD. 
For example dimensional regularisation is often 
used in perturbative studies and it breaks the 4- 
dimensional Lorentz invariance. Similarly lattice 
regularisation as proposed by Wilson Q breaks 
Lorentz invariance since a hypercubic lattice of 
space-time points is invoked. The key feature of 
Wilson's proposal is that gauge invariance is ex- 
actly retained. Then the approach to the con- 
tinuum limit (as the lattice spacing a is reduced 
to zero) can be shown to be well defined. Using 
Monte Carlo methods to explore lattice QCD, re- 
liable predictions can now be made for continuum 
quantities as I shall discuss. Moreover, in the con- 
tinuum limit, the Lorentz invariance is found to 
be fully restored. 

Lattice QCD needs as input the quark masses 
and an overall scale (conventionally given by 
Aqcd)- Then any Green function can be eval- 
uated by taking an average of suitable combina- 
tions of the lattice fields in the vacuum samples. 
This allows masses to be studied easily and ma- 
trix elements (particularly those of weak or elec- 



tromagnetic currents) can be extracted straight- 
forwardly. Unlike experiment, lattice QCD can 
vary the quark masses and can also explore dif- 
ferent boundary conditions and sources. This al- 
lows a wide range of studies which can be used to 
diagnose the health of phenomenological models 
as well as casting light on experimental data. 

It is worth recalling the weakness as well as 
the strength of the lattice approach to QCD. Any 
quantity which can be expressed as a vacuum ex- 
pectation value of fields can be extracted straight- 
forwardly in lattice studies. Thus masses and 
matrix elements of operators can be determined. 
What is not so easy is to explore hadronic de- 
cays. This is difficult because the lattice, using 
Euclidean time, has no concept of in and out 
states. About the only feasible strategies are to 
evaluate the mixing between states of the same 
energy - so giving some information on on-shcll 
hadronic decay amplitudes, or to make a model 
dependent analysis || of lattice results. 

One very special case is of considerable interest: 
this is quenched QCD where the sea-quark masses 
are taken as infinite. This suppresses quark loops 
in the vacuum completely, leaving just the full 
non-perturbative gluonic interactions. This glu- 
onic vacuum turns out to reproduce most of the 
salient features of QCD. It is also a very conve- 
nient approximation to use for comparison with 
phenomenological models. Quenched QCD is 
computationally rather easy to study and the pre- 
cise results allow the continuum limit to be ex- 
tracted reliably. Studies with sea quark effects 
included (known as dynamical fermion studies) 
are computationally much more demanding. I 
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discuss the current situation in this area in the 
last section. 

2. GLUEBALLS 

Glueballs are defined to be hadronic states 
made primarily from gluons. The full non- 
perturbative gluonic interaction is included in 
quenched QCD. In the quenched approximation, 
there is no mixing between such glueballs and 
quark - antiquark mesons. A study of the glue- 
ball spectrum in quenched QCD is thus of great 
value. This will allow experimental searches to be 
guided as well as providing calibration for models 
of glueballs. 

In principle, lattice QCD can study the meson 
spectrum as the sea quark mass is decreased to- 
wards experimental values. This will allow the 
unambiguous glueball states in the quenched ap- 
proximation to be tracked as the sea quark ef- 
fects are increased. It may indeed turn out that 
no meson in the physical spectrum is primarily 
a glueball - all states are mixtures of glue, qq, 
qqqq, etc. Studies conducted so far show no signif- 
icant change of the glueball spectrum as dynam- 
ical quark effects are added - but the sea quark 
masses used are still rather large |j| . 

In lattice studies, dimensionless ratios of quan- 
tities are obtained. To explore the glueball 
masses, it is appropriate to combine them with 
another very accurately measured quantity to 
have a dimensionless observable. Since the po- 
tential between static quarks is very accurately 
measured from the lattice (see the next section 
for more details), it is now conventional to use 
ro for this comparison. Here tq is implicitly de- 
fined by r 2 dV(r)/dr — 1.65 at r = ro where 
V(r) is the potential energy between static quarks 
which is easy to determine accurately on the lat- 
tice. In practice r may be related to the string 
tension a by ro^/a = 1.18. Taking the con- 
ventional value of y/a = 0.44 GeV then yields 
r = 2.68 GeV _1 =0.53 fm. Taking the potential 
V(r) from bb spectra suggests that ro w 0.5 fm 
which is a very similar estimate. 

Theoretical analysis indicates that for the 
quenched approximation the dimensionless ratio 
mro will differ from the continuum limit value by 
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Figure 1. The value of mass of the J PC = ++ 
and 2 ++ glueball states from refs[[|-[7| in units 
of ro where ro ~ 0.5 fm. The T2 and E repre- 
sentations are shown by octagons and diamonds 
respectively and their agreement indicates the 
restoration of rotational invariance for the 2 ++ 
state. The straight lines show fits describing the 
approach to the continuum limit as a — > 0. 



corrections of order a 2 . Thus in Fig. 1 the masses 
are plotted versus the lattice spacing a 2 for the 
J PC =0 ++ and 2 ++ glueballs. The straight lines 
then show the continuum limit obtained by ex- 
trapolating to a = 0. As can be seen, there is 
essentially no need for data at even smaller a- 
values to further fix the continuum value. The 
values shown correspond to m(0 ++ )r o = 4.33(5) 
and m(2 ++ )ro = 6.0(6). Since several lat- 
tice groups have measured these quan- 
tities, it is reassuring to see that the purely 
lattice observables are in excellent agreement. 
The publicised difference of quoted m(0 ++ ) from 
UKQCD §] and GF11 § comes entirely from re- 
lating quenched lattice measurements to values in 
GeV as I now discuss. 
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In the quenched approximation, different 
hadronic observables differ from experiment by 
factors of up to 10%. Thus using one quantity 
or another to set the scale, gives an overall sys- 
tematic error. Here I choose to set the scale by 
taking the conventional value of the string ten- 
sion, -y/er = 0.44 GeV, which then corresponds 
to Tq 1 = 373 MeV. An overall systematic error 
of 10% is then to be included to any extracted 
mass. This yields m(0++) = 1611(30)(160) MeV 
and m(2++) = 2232(220)(220) MeV where the 
second error is the systematic scale error. Note 
that these are glueball masses in the quenched ap- 
proximation - in the real world significant mixing 
with qq states could modify these values substan- 
tially. 

Recently a lattice approach using a large spa- 
tial lattice spacing with an improved action and 
a small time spacing has been used to study glue- 
ball masses. The results |^] arc that m(0 ++ )ro = 
3.98(15), m(2++)r = 5.85(2), m(l+-)r = 
7.21(2) and m'(2++)r = 8.11(4). There remains 
a small discrepancy with the result for the ++ 
glueball obtained above (m(0 ++ )r o = 4.33(5)) 
from lattice spacings much closer to the contin- 
uum limit. This has recently been explored in 
more detail ^ and modifications proposed to the 
technique which allow better control of the error 
in extrapolating to the continuum limit. 

I have focussed on the scalar and tensor glue- 
ball results because these are the lightest and 
best measured states in lattice studies. The glue- 
ball spectrum has been extracted for all J PC val- 
ues [|]]|]. Some of these results are shown in 
Fig. 2. One signal of great interest would be a 
glueball with J PC not allowed for qq - a spin- 
exotic glueball or oddball. These states are shown 
in Fig. 2 to be high lying: at least above 2m(0 ++ ). 
Thus they are likely to be in a region very difficult 
to access unambiguously by experiment. 

The only other candidate for a relatively light 
glueball is the pseudoscalar. Values quoted of 
m(0-+)r = 5.6(6), 7.1(1.1) and 5.3(6) from 
refs||H suggest an average of 6.0(1.0), close to 
the tensor glueball mass. This is confirmed by a 
study using asymmetric lattices which obtains 
a pseudoscalar glueball mass ratio to the tensor 
glueball of m(0- + )/m(2++) = 1.17(9). 



Within the quenched approximation, the glue- 
ball states are unmixed with qq, qqqq, etc. Fur- 
thermore, the qq states are 'magically mixed' in 
the quenched approximation. Once quark loops 
are allowed in the vacuum, for the favour-singlet 
states of any given J PC , there will be mixing be- 
tween the ss state, the uu + dd state and the 
glueball. For the ++ case, as a guide to this 
mixing, the quenched mass spectrum of qq states 
has been determined on a lattice and the scalar 
mesons are found to lie somewhat lighter than 
the tensor states ]io|] . These 2 ++ mesons are ex- 
perimentally approximately 'magically mixed' so 
will be quite close to the quenched mass deter- 
mination. This suggests that the quenched scalar 
masses from the lattice are at around 1.2 GeV 
and 1.5 GeV. An independent study QJl^l sug- 
gests that the scalar ss state is about 200 MeV 
lighter than the glueball which is a broadly com- 
patible conclusion. Thus the glueball, at around 
1.6 GeV, lies heavier than the qq scalar states. 

A first attempt Jl^] has been made to measure, 
on a quenched lattice, the mixing matrix element 
between the glueball and the scalar qq states. The 
mixing matrix element between the ss meson and 
the glueball is extracted as E(s) = 43(31) MeV. 
The mixing matrix elements with the glueball are 
found to be similar for light and ss quarks (actu- 
ally E(u, d)/E(s) ps 1.2). With this lattice input, 
a mixing model |l^,|l^| allows the content of the 
experimental scalar states at 1390, 1500 and 1710 
MeV to be assigned. Other, more phenomenolog- 
ical mixing models have also been proposed Jl4[ . 
In these mixing scenarios none of the experimen- 
tal states is predominantly an unmixed glueball. 

As well as this mixing of the glueball with qq 
states, there will be mixing with qqqq states which 
will be responsible for the hadronic decays. Be- 
cause of the limitations of lattice techniques, it is 
only possible to estimate these decay matrix ele- 
ments when the initial and final states have the 
same energy and this has been applied to the de- 
cay matrix element between a ++ glueball and 
two pseudoscalar qq mesons. This then allows an 
estimate of the width which would arise when go- 
ing beyond the quenched approximation. A first 
attempt to study this jl^] yields estimated widths 
of order 100 MeV. Even though this lattice cal- 
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Figure 2. The mass of the glueball states with quantum numbers J PC from ref[^|. The scale is set by 
yfd w 440 MeV which yields the right hand scale in MeV. The solid circles represent mass determinations 
whereas the open circles are upper limits. Note that excited states in several cases have been determined, 
such as m'(0++)/V^ = 6.3(4) from refg. 



culation is very exploratory, it does indicate that 
very wide widths to two pseudoscalars are not ex- 
pected. A more realistic study would involve tak- 
ing account of mixing with the nn and ss scalar 
mesons as well as the decay channels. 

3. HEAVY QUARK INTERACTIONS 

In the limit tuq — > oo, the heavy quark effective 
theory describes a universal behaviour. For finite 
rriQ , corrections of order 1 /tuq are expected. The 
simplest way to study the heavy quark limit on 
a lattice is to use static quarks. The potential 
energy V(R) between a static quark and anti- 
quark at separation R is readily obtained. Then 
for heavy quarks, one may solve for the spectra 



in this potential using the Schrodinger equation 
in the adiabatic approximation. This is expected 
to be a good approximation for determining the 
spectrum of bb mesons. The quenched lattice po- 
tential is well measured and is found to have a 
form parametrised by 



V(R) = V --+aR 
it 



(1) 



where e is the coefficient of the Coulomb term and 
a is the string tension. This expression shows 
that the potential continues to increase as R is 
increased - this is confinement. 

A comparison from ref[fl6| of the spectrum in 
the quenched lattice potential with the T states 
is shown in Fig. 3. The lattice result is quali- 
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Figure 3. Potentials V(R) in GeV between static 
quarks at separation R (in units of tq w 0.5 fm) 
for the ground state (square and *) and for the E u 
symmetry which corresponds to the first excited 
state of the gluonic flux (octagons and diamonds). 
Results from the quenched calculations of refl]L6| 
are shown by symbols corresponding to different 
lattice spacings. For the ground state potential 
the continuous curve is an interpolation of the 
lattice data while the dotted curve with enhanced 
Coulomb term fits the spectrum and yields the 
masses shown. The lightest hybrid level in the 
excited gluonic potential is also shown. 



tatively similar to the experimental T spectrum. 
The main difference is that the Coulombic part 
(e) is effectively too small (0.28 rather than 0.48). 
This produces jl6) a ratio of mass differences 
(IP—1S)/(2S— IS) of 0.71 to be compared with 
the experimental ratio of 0.78. This difference in 
e is understandable as a consequence of the one 
gluon exchange (Coulombic) force at short dis- 
tances which would be increased by 33/ (33— 2 A^) 
in perturbation theory in full QCD compared to 
quenched QCD. We will return to discuss this. 



The situation of a static quark and antiquark 
is a very clear case in which to discuss hybrid 
mesons which have excited gluonic contributions. 
A discussion of the colour representation of the 
quark and antiquark is not useful since they are 
at different space positions and the combined 
colour is not gauge invariant. A better criterion 
is to focus on the spatial symmetry of the gluonic 
flux. As well as the symmetric ground state of 
the colour flux between two static quarks, there 
will be excited states with different symmetries. 
These were studied on a lattice [|l7| and the con- 
clusion was that the E u symmetry (corresponding 
to flux states from an operator which is the differ- 
ence of U-shaped paths from quark to antiquark 
of the form PI — U) was the lowest lying gluonic 
excitation. Results for this potential are shown 
in Fig. 3. 

This gluonic excitation in the E u representa- 
tion (or n u in the continuum) corresponds to 
a component of angular momentum of one unit 
along the quark antiquark axis. Then one can 
solve for the spectrum of hybrid mesons using 
the Schrodinger equation in the adiabatic ap- 
proximation. The spatial wave function neces- 
sarily has non zero angular momentum and the 
lightest states correspond to L PC = 1 H and 
1 K Combining with the quark and antiquark 
spins then yields Jl7j a set of 8 degenerate hybrid 
states with J PC = 1~, 0"+, 1"+, 2 -+ and 
1 ++ , H , 1 H , 2 H respectively. These con- 
tain the spin-exotic states with J PC = 1 H 
and 2 H which will be of special interest. 

Since the lattice calculation of the ground state 
and hybrid masses allows a direct prediction for 
their difference, the result for this 8-fold degener- 
ate hybrid level is illustrated in Fig. 3 and corre- 
sponds Q to masses of 10.81(25) GeV for bb and 
4.19(15) GeV for cc when using the Schrodinger 
equation in the Born-Oppcnhcimcr approxima- 
tion to determine the spectrum. Here the errors 
take into account the uncertainty in setting the 
ground state mass using the quenched potential 
as discussed above. Recently a different lattice 



technique 18 has been used to explore in detail 
the excited gluonic levels in the quenched approx- 
imation. The results above are confirmed and 
the value quoted for the lightest hybrid meson 
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is (mnybrid — TOr^o = 3.166(3) for 66 where the 
error does not take into account unquenching 

The quenched lattice results, after adjusting to 
take account of the measured 66 spectrum, sug- 
gest that the lightest hybrid mesons lie above the 
open BB threshold by about 270 MeV. This can 
also be studied by comparing directly the lattice 
hybrid masses with twice the quenched lattice 
masses for the B meson Q. Using quenched 
results from the smallest lattice spacing (/3 = 
6.2) available with clover-improved fermions pQ | 
yields m Hy brid ~ 2m s w 140(80) MeV. This es- 
timate is somewhat smaller than that obtained 
above. In both cases, however, the hybrid lev- 
els lie above the open threshold and are likely 
to be relatively wide resonances. Another conse- 
quence |ll| is that the very flat potential implies 
a very extended wavefunction. This is illustrated, 
for example, in ref jl8| which shows that the hy- 
brid wavefunction remains significant out to a ra- 
dius of 1 fm. This has the implication that the 
wavefunction at the origin will be small, so hy- 
brid vector states will be weakly produced from 

It would be useful to explore the splitting 
among the 8 degenerate J PC values obtained. 
This could come from different excitation energies 
in the L PC = 1 H (magnetic) and 1 h (pseudo- 
electric) gluonic excitations, spin-orbit terms, as 
well as mixing between hybrid states and QQ 
mesons with non-exotic spin. One way to study 
this on a lattice is to use the NRQCD formu- 
lation which describes non-static heavy quarks 
which propagate non-relativistically. Prelimi- 
nary results for hybrid excitations from several 
groups [ pTp^Jl^ l give at present similar results 
to those with the Born-Oppenheimer approxi- 
mation in the static approximation as described 
above. There is some evidence that the mag- 
netic excitations (which include the 1 ^ spin ex- 
otic) are lighter than the pseudo-electric ones 
(which include the 0" 1 spin exotic). For exam- 
ple ref@ obtains (m(l+ - ) - rn T )r = 3.29(5) 
and (m(0 '") — mx)ro = 3.51(8), to be com- 
pared with the static approximation, discussed 
above, which has degenerate levels with a value of 
3.166(3). It has to be remembered that these are 
all quenched determinations, so there are system- 



atic errors of order 10% from unquenching. This 
implies that the NRQCD estimate for the light- 
est spin-exotic hybrid would be the 1 H state at 
a mass of 10.69(13) GeV. 

4. LIGHT QUARK INTERACTIONS 

Unlike very heavy quarks, light quark propaga- 
tion in the gluonic vacuum sample is very compu- 
tationally intensive - involving inversion of huge 
(10 9 x 10 9 ) sparse matrices. Current computer 
power is sufficient to study light quark physics 
thoroughly in the quenched approximation. The 
state of the art |23| is the Japanese CP-PACS 
Collaboration who are able to study a range of 
large lattices (up to about 64 4 ) with a range 
of light quark masses. Qualitatively the meson 
and baryon spectrum of states made of light and 
strange quarks is reproduced with discrepancies 
of order 10% in the quenched approximation. 

Here I will focus on hybrid mesons made from 
light quarks. There will be no mixing with qq 
mesons for spin-exotic hybrid mesons and these 
are of special interest. The first study of this area 
was by the UKQCD Collaboration |10| who used 
operators motivated by the heavy quark studies 
referred to above. Using non-local operators, they 
studied all 8 J PC values coming from L PC = 1 H 
and 1 h excitations. The resulting mass spec- 
trum is shown in Fig. 4 where the J PC — 1 
state is seen to be the lightest spin-exotic state 
with a statistical significance of 1 standard devi- 
ation. The statistical error on the mass of this 
lightest spin-exotic meson is 7% but, to take ac- 
count of systematic errors from the lattice de- 
termination, a mass of 2000(200) MeV is quoted 
for this hybrid meson with ss light quarks. Al- 
though not directly measured, the corresponding 
light quark hybrid meson would be expected to 
be around 120 MeV lighter. In view of the small 
statistical error, it seems unlikely that the 1 h 
meson in the quenched approximation could lie 
as light as 1.4 GeV where there are experimental 
indications for such a state |2~jf| . Recently exper- 
imental evidence has also been presented for 
another 1 h meson at 1.6 GeV. 

One feature clearly seen in Fig. 4 is that non 
spin-exotic mesons created by hybrid meson oper- 
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Figure 4. The masses in GeV of states of J PC 
built from hybrid operators with strange quarks, 
spin-exotic (*) and non-exotic (squares). The 
dot-dashed lines are the mass values found for 
ss operators. Results from refflQ]. 



ators have masses which are very similar to those 
found when the states are created by qq oper- 
ators. This suggests that there is quite strong 
coupling between hybrid and qq mesons even in 
the quenched approximation. This would imply 
that the 7r(1800) is unlikely to be a pure hybrid, 
for example. 

A second lattice group has also evaluated 
hybrid meson spectra from light quarks from 
quenched lattices. They obtain [£6| masses of the 
1 h state with statistical and various systematic 
errors of 1970(90)(300) MeV, 2170(80)(100)(100) 
MeV and 4390(80) (200) MeV for nn, ss and cc 
quarks respectively. For the 0" 1 spin-exotic state 
they have a noisier signal but evidence that it 
is heavier. They also explore mixing matrix el- 
ements between spin-exotic hybrid states and 4 
quark operators. 

Recently a first attempt has been made [E7| to 



determine the hybrid meson spectrum using full 
QCD. The sea-quarks used have several different 
masses and an extrapolation is made to the limit 
of physical sea-quark masses, yielding a mass of 
1.9(2) GeV for the lightest spin-exotic hybrid me- 
son, which they again find to be the 1 In 
principle this calculation should take account of 
sea-quark effects such as the mixing between such 
a hybrid meson and qqqq states such as t]tt. As il- 
lustrated in Fig. |[ the calculations are performed 
for quite heavy sea quarks (the lightest being ap- 
proximately the strange quark mass) and then a 
linear extrapolation is made. It is quite possible, 
however, that such mixing effects turn on non- 
linear ly as the sea-quark masses are reduced. The 
systematic error from this possibility is difficult to 
quantify. 

The three independent lattice calculations of 
the light hybrid spectrum are in good agreement 
with each other. They imply that the natu- 
ral energy range for spin-exotic hybrid mesons is 



around 1.9 GeV. The J 
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state is found 



to be lightest. It is not easy to reconcile these 
lattice results with experimental indications for 
resonances at 1.4 GeV and 1.6 GeV. Mixing with 
qqqq states such as rjiT is not included for realistic 
quark masses in the lattice calculations. This can 
be interpreted, dependent on one's viewpoint, as 
either that the lattice calculations are incomplete 
or as an indication that the experimental states 
may have an important meson-meson component 
in them. 

5. TOWARDS FULL QCD 

Algorithms exist which allow lattice simulation 
of full QCD with sea quarks of mass m sea . This 
study needs lots of CPU power since the sea quark 
loops in the vacuum are represented effectively 
as a long range interaction between the gluonic 
degrees of freedom. Most studies to date have 
been exploratory with sea quark masses above the 
strange quark mass. In this regime, very little 
change from the quenched approximation is seen 
in physical predictions from the lattice. 

One area where specific changes are expected 
is for the potential between heavy quarks. At 
small separation R, the Coulombic term is ex- 
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Figure 5. The extrapolation in sea-quark mass 
(in lattice units) from ref[^7j for the 1 h hybrid 
meson. The dotted vertical line corresponds to 
light sea quarks. 



pected to increase in strength and indeed evi- 
dence of this has been reported 28 - 30| . At larger 
R, signs of string breaking are expected since a 
light quark antiquark pair can be produced from 
the vacuum to yield two Qq mesons with energy 
independent of R at large separation. This has 
been explored by the SESAM, UKQCD and CP- 
PACS collaborations and little sign of the effect 
is seen gUg. 

It has been very difficult to see unambiguous 
signs of sea quark effects in the hadron spec- 
trum. The CP-PACS collaboration |l) have re- 
cently reported that they see evidence that, for 
two flavours of light quark, results are closer to 
experiment than from quenched studies. Even 
their study, however, is restricted to sea-quark 
masses which are not appreciably lighter than the 
strange quark mass. Now, it is possible that such 
sea-quark effects turn on non-linearly as the sea 
quark mass is reduced. As an example, in cur- 



rent studies the 7r + 7r P-wave is heavier than the 
p meson so the p cannot decay. Another example 
is that glueball mixing with qq scalar mesons may 
depend quite sensitively on their masses, so as the 
sea quark mass is reduced, the scalar meson mass 
will move relative to the glueball mass so causing 
a non-linear change in the mixing strength. 

Further work is needed to reduce the sea quark 
mass and to increase the lattice size. Dedicated 
computing power of several hundreds of Gfiops is 
available to lattice collaborations and, with future 
plans for Tflops machines, progress in this area 
should now be possible. 



6. SUMMARY 

Quenched lattice QCD is well understood 
and accurate predictions in the continuum limit 
are increasingly becoming available. Glue- 
ball masses of m(0++) = 1611(30)(160) MeV; 
m(2++) = 2232(220)(220) MeV and m(0"+) = 
2611(370)(260) MeV are predicted where the sec- 
ond error is an overall scale error. The quenched 
approximation also gives information on quark- 
antiquark scalar mesons and their mixing with 
glucballs. This determination of the mixing in the 
quenched approximation is currently more power- 
ful than attempting to determine the mixed spec- 
trum directly in full QCD because the sea-quark 
masses realisable are rather heavy. There is also 
some lattice information on the hadronic decay 
amplitudes of glueballs. 

For hybrid mesons, there will be no mixing with 
qq for spin-exotic states and these are the most 
useful predictions. The J PC = 1 h state is ex- 
pected at 10.81(25) GeV for b quarks; 4.19(15) 
GeV for c quarks, 2.0(2) GeV for s quarks and 
1.9(2) GeV for u, d quarks. Mixing of spin-exotic 
hybrids with qqqq or equivalently with meson- 
meson is allowed and will modify the predictions 
from the quenched approximation. 

Much activity is currently underway to explore 
the effects of sea quarks of ever decreasing mass. 
Future teraflops computing facilities will be es- 
sential to obtain quantitative results for hadronic 
spectroscopy. 
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